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ABSTRACT 

The infrared detection of a z > 7 quasar (Mortlock et al., 201 1) has opened up a window to 
, directly probe the inter-galactic medium (IGM) during the epoch of reionization. It is antic- 

ipated that future observations will yield more quasars extending to higher redshifts. In this 
paper we theoretically consider the possibility of detecting the ionized bubble around a z = 8 
quasar using targeted redshifted 21 -cm observations with the GMRT. The apparent shape and 



o 



O " size of the ionized bubble, as seen by a distant observer, depends on the parameters Nph,s/C, 

^ xh I /C and tq where Nphs and tq are respectively the ionizing photon emission rate and 



C/3 



in 



age of the quasar, and xh i and C are respectively the neutral fraction and clumping factor of 
the IGM. The 21 -cm detection of an ionized bubble, thus, holds the promise of allowing us to 
probe the quasar and IGM properties at z = 8. 
'nJ" ■ In the current work we have analytically calculated the apparent shape and size of a 

^ . quasar's ionized bubble assuming an uniform IGM and ignoring other ionizing sources be- 

sides the quasar, and used this as a template for matched filter bubble search with the GMRT 
visibility data. We have assumed that Nphs is known from the observed infrared spectrum, 
\^ • and C = 30 from theoretical considerations, which gives us the two free parameters x-a i 

and Tq for bubble detection. Considering 1, 000 hr of observation, we find that there is a rea- 
sonably large region of parameter space bounded within {x-a i , (rg/lO^yr)) = (1.0, 0.5) 
and (0.2, 7.0) where a 3a detection is possible if (iVp/is/10^'' sec^^) = 3. The available 
region increases if Nphs is larger, whereas we need xh i > 0.4 and (rg/lO^yr) > 2.0 if 
(-/Vp/is/10^'' sec~^) = 1.3. Considering parameter estimation, we find that in many cases it 
will be possible to quite accurately constrain tq and place a lower limit on xh i with 1 , 000 
hr of observation, particularly if the bubble is in the early stage of growth and we have a very 
luminous quasar or a high neutral fraction. Deeper follow up observations (4, 000 and 9, 000 
hr) can be used to further tighten the constraints on tq and .th i ■ We find that the estimated 
I is affected by uncertainty in the assumed value of C. The quasar's age tq however is 
robust and is unaffected by the uncertainty in C. 

The presence of other ionizing sources and inhomogeneities in the IGM distort the shape 
and size of the quasar's ionized bubble. This is a potential impediment for bubble detection 
and parameter estimation. We have used the semi-numerical technique to simulate the ap- 
parent shape and size of quasar ionized bubbles incorporating these effects. If we consider a 
9, 000 hr observation with the GMRT we find that the estimated parameters tq and xh i are 
expected to be within the statistical uncertainties. 

Key words: methods: data analysis - cosmology: theory: - diffuse radiation 



1 INTRODUCTION 



The epoch of reionization is one of the least known periods in 
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2011). The intergalactic medium (IGM) during this period is char- 
acterised by bubbles of ionized hydrogen (H II ), centred around 
luminous sources. Stars forming within early galaxies are believed 
to be the major sources of ionizing photons during this era (for re- 
views see Choudhury & Ferrara 2006; Choudhury 2009). Quasars 
are expected to be very rare during this epoch, but they are capa- 
ble of generating larger H II bubbles around them with respect to 
their stellar counterpart. Detection of these H II bubbles will di- 
rectly probe the state of the local IGM around the ionizing sources. 
It will also constrain the properties of the quasar, such as its lumi- 
nosity and age. 

The highest redshift at which a quasar has been detected till 
date is 2 = 7.085 (Mortlock et al., 2011). The Ly-a absorption 
spectrum of this object reveals a highly ionized near zone of radius 
2.1 ± 0.1 Mpc (physical). The observed near zone size is unexpect- 
edly small when compared with the other high redshift quasars ob- 
served in the redshift range 6.0 < z < 6.4 (Fan et al., 2003; Willott 
et al., 2007, 2010). The ionizing photon emission rate estimated for 
this quasar is 1.3 x lO^^sec^^ assuming a power law for the un- 
absorbed continuum emission blueward of Ly-a line. The quasar's 
age (tq) and local volume averaged neutral fraction ((sh i )v) can 
be estimated from the size of the quasar's near zone, however the 
estimated tq and (a;H i )v are strongly correlated. For the same ob- 
servation Bolton et al. (201 1) have constrained the tq and (xh i )v 
using simulated Ly-a absorption spectra. This study shows that 
there are several combinations of the tq and (sh i )v all of which 
can reproduce the observed data. They find that the spectrum ob- 
served by Mortlock et al. (2011) is consistent with tq ~ 10^ yr 
with either (ihi)^ ~ 10^'* — 10^'^ provided there is a prox- 
imate damped Ly-a absorber (DLA) or (xhi)v' ~ 0.1 without 
a proximate DLA. The same observation is also consistent with 
TQ ~ 10^ yr and either (xh i )v < 10^* with a proximate DLA 
or a fully neutral IGM (xh i )v ~ 1 without a proximate DLA. 
It is quite clear that observations of Ly-a absorption spectra are 
limited in their ability to determine the age of the quasar and the 
local neutral fraction. This limitation essentially arises from two 
reasons. The first being that the IGM becomes nearly opaque at a 
very low neutral fraction (xh i)v ~ 10"'*. It is not possible to 
distinguish between a neutral fraction of 10"* from a completely 
neutral medium, and the actual ionized bubble around the quasar 
may extend far beyond the region inferred from the Ly-a absorp- 
tion spectra. The second limitation arises from the fact that the Ly- 
a absorption spectrum is a pencil beam observation along a line of 
sight (LoS) to the quasar. The presence of a proximate DLA can 
completely change the interpretation of the spectrum. 

The redshifted 21-cm emission from neutral hydrogen in the 
epoch of reionization is believed to be a powerful tool to detect 
H II bubbles around quasars. The intensity is directly proportional 
to the H I density, and it is in principle possible to probe the entire 
ionization profile of the H II bubble. It may be possible to over- 
come the limitations of Ly-a absorption spectra using redshifted 
21-cm observations, and place better constraints on the quasar pa- 
rameters and the state of the IGM. This is particularly motivated 
by the fact that the presently functioning GMRT* (Swarup et al., 
1991) and LOFAR^ and the upcoming MWA'^ and 2 1 CM A ■* are 
all sensitive to the H I signal from the epoch or reionization. How- 
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ever, redshifted H I observations have their own limitations in that 
the H I signal is extremely small relative to the sensitivities of the 
present and upcoming telescopes. Further, the signal will be buried 
deep in foregrounds which are a few orders of magnitude larger 
than the signal (Ali, Bharadwaj & Chengalur, 2008; Bernardi et al., 
2009). 

Datta, Bharadwaj & Choudhury (2007) (hitherto Paper I) have 
proposed a matched filter technique to detect ionized bubbles in 
radio-interferometric observations of redshifted 21-cm emission. 
The matched filter optimally combines the entire signal of an ion- 
ized bubble while minimizing the noise and the foreground contri- 
butions. The redshift z ~ 8 is most optimal for bubble detection 
(Datta, Bharadwaj & Choudhury, 2009), and at this redshift a 3cr 
detection is possible for bubbles of comoving radius larger than 
24 and 28 Mpc with 1000 hr of observation using the GMRT and 
MWA respectively, assuming that the neutral hydrogen fraction is 
~ 0.6 outside the bubble. This technique, however, is limited by 
the density fluctuations in H I outside the bubble, and analytical es- 
timates (Paper I) and simulations (Datta et al., 2008) show that it is 
not possible to detect bubbles of comoving radius < 6 and < 12 
Mpc using the GMRT and MWA respectively, however large be the 
observation time. 

The matched filter technique mentioned above assumes the 
H II bubble to be spherical. However, a growing spherical bubble 
will appear anisotropic along the LoS to a present day observer due 
to the finite light travel time (FLTT) (Wyithe & Loeb, 2004; Yu, 
2005; Wyithe, Loeb & Barnes, 2005; Sethi & Haiman, 2008) and 
also due to the evolution of global ionization fraction (Geil et al., 
2008). In an earlier work (Majumdar et al. (2011); hitherto Paper 
II), we have analytically quantified this anisotropy and studied the 
possibility of detecting such a bubble in a targeted search around 
a known quasar. We find that the bubble appears elongated along 
the LoS during the early stages of its growth, whereas it appears 
compressed in the later stages when the growth slows down. In ad- 
dition to this, the apparent centre of the bubble also shifts towards 
the observer. We find that a spherical filter is adequate for bubble 
detection even when the apparent anisotropy of the bubble's shape 
is taken into account, the centre of the best matched filter will how- 
ever be shifted relative to the quasar. We also propose that the mea- 
sured shift and the radius of the best matched filter can together be 
used to constrain the age and luminosity of the quasar. 

In this work we consider the detection of an H II bubble in a 
targeted 21-cm search centred on a known quasar at 2 = 8 and 
we investigate the possibility of using such an observation to con- 
strain the properties of the quasar and the IGM. For this purpose 
we improve upon the spherical filter by calculating the apparent, 
anisotropic shape of the quasar bubble and using this as a template 
for the filter. We expect this to give a better match to the bubble that 
is actually present in the data. Further, we also expect improved 
parameter estimation using the improved filter. 

Our earlier work (Paper II) was entirely based on analytic es- 
timates which do not take into account the presence of ionizing 
sources other than the quasar and the inhomogeneities in the IGM. 
We have overcome this limitation here by using the semi-numerical 
technique (Choudhury, Haehnelt & Regan, 2009) to simulate the 
ionization field. 

This paper is arranged as follows - In Section 2 we briefly dis- 
cuss the model for bubble growth around a quasar, and the matched 
filter technique for detecting such a bubble. We also present the im- 
proved filter based on the calculated apparent shape of the bubble. 
In Section 3 we use analytic estimates, based on the bubble growth 
equation, to study the SNR for bubble detection. We use this to 
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Figure 1. This shows the apparent shape of the H II bubble around a quasar 
Q as seen by a distant observer located in the direction indicated by the 
arrows. The apparent anisotropy of the H II bubble arises from the fact that 
photons received from different parts of the bubble had been emitted at 
different stages of of its growth. The apparent centre of the bubble is shifted 
by QC from the quasar Q to the point C. The bubble's apparent radius is 
EC = CB = Ri] along the LoS and CD = ij^ perpendicular to LoS. 



determine the parameter range where a 3ct detection is possible in 
1, 000 hr of observation with the GMRT. We then consider param- 
eter estimation, and explore the kind of limits that can be placed on 
the quasar and IGM properties using the matched filter technique. 
Section 4 briefly describes how we have simulated the apparent 
shape of quasar bubbles, and we present our results and summarize 
our findings in Section 5. 

Our entire analysis is restricted to the redshift 2 = 8. We have 
used the values from the WMAP 7 year data (Komatsu et al., 2011 ; 
Jarosik et al., 2011) h = 0.705, Vim = 0.2726, = 0.726, 
Q,bh? — 0.0223 for the cosmological parameters. 



2 BUBBLE DETECTION 
2.1 Model for Bubble Growth 

The growth of a spherical H II region around a quasar which is 
isotropically emitting ionizing photons at a rate NpUa is governed 
by the equation (Shapiro & Giroux, 1987; White, Becker, Fan & 
Strauss, 2003; Wyithe & Loeb, 2004; Wyithe, Loeb & Barnes, 
2005; Yu, 2005; Sethi & Haiman, 2008), 
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[eq. (1) of Paper II and eq (7) of Yu 2005]. It is assumed that the 
quasar is triggered at a cosmic time t;, and r = t — ti denotes the 
quasar's age at any later time t. The variable r (r ) denotes the radius 
of the spherical ionizing front at the instant when a photon that was 
emitted from the quasar at r catches up with the ionizing front. The 
term as(= 2.6 x 10~^^ cm'^ s^^) is the recombination coefficient 
to excited levels of hydrogen at T = lO** K, {nn} and {nn i ) are 
the average hydrogen and neutral hydrogen densities respectively, 
xh I is the neutral fraction of the IGM and C = {n%i , ) / {uhY 
is the clumping factor which quantifies the effective dumpiness 
of the hydrogen inside the bubble. Eq. (1) tells us that the growth 
of the H II bubble is determined by the ionizing photon emission 
rate after accounting for the photons required to compensate for 
the recombinations inside the existing ionized region. 



For a constant NpUs 
(eq. (1)) takes the form 

r(r) = rs 



the solution to the above growth model 



exp 



[eq. (2) of Paper II and eq. (8) of Yu 2005] where Tr 
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I (C (nn) oibY and rs = [^NphsTvEc/ (jiiix-a i {nH))j 
The apparent shape of the growing bubble that is perceived by a 
present day observer will be distorted due to the fact that photons 
received by the observer from different parts of the bubble had been 
emitted at different stages of the bubble's growth (see Figure 1). To 
visualize this apparent shape one need the relation between r and 
the angle between the LoS and the point A (Figure 1) under con- 
sideration on the ionization front. Consider a quasar observed at 
an age tq. The light travel time starting from the quasar at r to 
the point A and then to the observer is [r(r)/c](l — cos.(f)) more 
compared to the photon that was emitted from the quasar at tq and 
travels straight to the observer. This gives 



c 



(3) 



[eq. (7) of Paper II and eq. (3) of Yu 2005]. We use eq. (2) and 
eq. (3) to determine r as a function of 0. This gives the quasar's 
apparent shape shown in Figure 1 . 

The apparent shape of a quasar's H II bubble has been studied 
in detail in Paper II. We find that the FLTT has two effects, 1. the 
bubble appears to be anisotropic along the LoS and 2. the bubble's 
apparent centre shifts along the LoS away from the quasar towards 
the observer. We now briefly discuss how we parametrize these two 
effects. Referring to Fig 1, we use R±_ = CD to quantify the over- 
all comoving size of the bubble. The bubble centre C is located 
mid way between EB. We use the dimensionless shift parameter s 
defined as. 



QC 



(4) 



to quantify the shift in the bubble's centre as a fraction of its radius. 
We find (see Paper II for details) that a can be greater than 1 during 
the early stages of the bubble's growth and it approaches in the 
later stages of evolution. 

We use the dimensionless anisotropy parameter rj defined as. 



CB 



1 , 



(5) 



to quantify the anisotropy in the apparent shape of the bubble. A 
value 77 > indicates that the bubble is elongated along the LoS 
and a value 77 < 0, indicates that it is compressed along the LoS. 
Our earlier work (Paper II) shows that rj has values in the range 0.1 
to 0.5 during the early stages of growth when the bubble appears 
elongated along the LoS. We also find that ry has values in the range 
0.0 to —0.2 in the later stages of evolution when the bubble may 
appear compressed along the LoS. 



2.2 Matched filter bubble detection 

The basic observable quantity in the radio interferometry is the vis- 
ibility V{U , v) which is related to the specific intensity pattern on 
the sky (6) as 



V{U,y) 



(6) 
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The baseline U — d/X denotes the antenna separation d projected 
in the plane perpendicular to the LoS in units of the observing 
wavelength \, 9 is a two dimensional vector in the sky plane with 
the origin at the centre of the FoV (the phase centre) and A{9) 
is the beam pattern of a single antenna. For the GMRT ^(6^) can 
be well approximated by a Gaussian A{6) = e~® where 
6*0 « 0.6 Sfwhm and we use 2.28° for 6*0 at 157.77 MHz i.e. 
z = 8. We consider a situation where the observation is spanned 
across several frequency channels around the central frequency of 
157.77 MHz. 

The visibility recorded in a radio-interferometric observation 
is actually a combination of several contributions, 

V{U, v) = S0, v) + HFijj, v) + v) + F(tj, v) , (7) 

where the first term S{U , v) is the signal from the ionized region 
that is actually present in the observational data, HF{U, v) is the 
contribution from fluctuations in the H I distribution outside the 
ionized bubble, N{U, v) is the system noise which is inherent to 
the measurement and F{U , u) is the contribution from astrophysi- 
cal foregrounds. 

The signal S{U , u) from an ionized region will appear as 
a decrement with respect to the background 21-cm radiation. A 
spherical H II bubble embedded in an uniform IGM with a neu- 
tral fraction xh i is parametrized by its comoving radius Ri,, the 
redshift Zc and the angular position 6c corresponding to the centre 
of the bubble. The quantity is the comoving distance to the red- 
shift where the H I emission received al i/ = 1420 MHz/(l + z) 
originated, and rl = dri^/du. A plane perpendicular to LoS and 
passing through the centre of the bubble at a comoving distance 
will cut a disk of comoving radius i?^ = Rt a/I — {Au / Aui,)^ 
from the bubble, where Av = Vc — v and Avh — Rt/r'^^, and 
6v = Rv/ti, is the angular extent corresponding to Rv. For a bub- 
ble located at the centre of the FoV we have. 
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where L = 2.5 x 10^ (£^) ((-I) (_|^), e(x) is the 

Heaviside step function, and Ji (x) is the Bessel function of first 
order. The expected signal has a peak value (Paper I) of ^(O, i^) — 
1.12 mjy for a bubble with Ri, = 40 Mpc, and S{0, v) <x Rl. 

The expected H I signal is very weak, so much so that the con- 
tribution HF(U, v) from H I fluctuations outside the bubble may, 
in some cases, exceed the expected signal. We consider the system 
noise N{U, u) in each baseline and frequency channel to be an in- 
dependent Gaussian random variable with zero mean and r.m.s. 
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where the constant is different for different interferometric ar- 
rays, Ai/ is the channel width and At is the correlator integration 
time. Following Paper I, we use C = 1.03 Jy for the GMRT. We 

note that the r.m.s. noise 



(N^) will actually vary with U, we 
have ignored this baseline dependence in order to keep the analysis 
simple. Further we do not expect this to be a very significant effect 
compared to the various other uncertainties arising from the lack 
of knowledge about the IGM and the quasar properties. The astro- 
physical foregrounds F{U, v) are expected to be several orders of 
magnitude larger than the H I signal. But they are predicted to have 



a featureless, continuum spectra whereas the signal is expected to 
have a dip around Pc (central frequency of the target H II bubble). 
This difference holds the promise of allowing us to separate the 
signal from the foregrounds. 

The visibility based matched filter technique, proposed in Pa- 
per I, optimally combines the signal from an ionized region while 
minimizing the other contributions. In this technique we use the 
signal expected from a spherical H II bubble centered at redshift Zb 
and comoving radius R^ as a template, and search for the presence 
of this signal in the data using the estimator, 



(10) 



Here Sf{U, u) is the filter and V{U , u) are the measured visibili- 
ties. The filter is defined as. 
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where S{U , v) is the signal expected from the bubble that we are 
trying to detect. The filter eliminates any frequency independent 
component of the foreground from the frequency range Vc + B' /2 
to Uc — B' 12. All the contributions to V, except the signal SiJJ , v), 
are assumed to be random variables of zero mean, uncorrelated to 
the filter whereby the estimator has expectation value. 



(12) 



The other terms in eq. (7) contribute only to the variance, 

{{AEf) = ((AS)')hf + ((A^)')n + ((A£;)")f . (13) 
The signal to noise ratio for the estimator is 



SNR: 



{{AEr) 



(14) 



Matched filter bubble detection is carried out by analyzing the SNR 
for different values of the filter parameters Zf and Rf. We ex- 
pect the SNR to peak when the filter parameters Zf and Rj ex- 
actly match the parameters zt and Rb of the bubble that is actually 
present in the data. We have a statistically significant (3a) bubble 
detection if the peak SNR > 3. Note that the two angular coordi- 
nates 9x and 9y are already known in a targeted search around a 
quasar. In general, 9x and 9y have also to be treated as free param- 
eters in a blind search. 

The filter Sf{U, v) is defined in such a way [eq. (11)] that 
it is insensitive to the presence of a smooth frequency indepen- 
dent foreground component. Under the assumed foreground model 
(Paper I), the residual foreground contribution {{AE)'^)y is pre- 
dicted to be much smaller than the other contributions to {{AE)"^) 
and we do not consider it here. In the absence of patchy reioniza- 
tion outside the bubble, the H I fluctuations trace the dark matter 
fluctuations. This imposes a lower bound Rt > 12 Mpc for the 
GMRT on the smallest H II bubble that can be detected (Datta et 
al., 2008). Considering patchy reionization outside the bubble, we 
expect galaxies to produce ionized regions with a typical radius of 
6 Mpc or smaller. The quasar bubbles that we can detect using the 
GMRT (> 20 Mpc) are much larger than these ionized patches. 
Earlier work (Datta et al., 2008) shows that the H I fluctuations do 
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Figure 2. This shows the H II bubble around a quasar with scc^^) = 3 at different stages of growth ((tq /lO'' yr) = 0.5, 2.0, 4.0). The three 

different shapes in each panel con'espond to different IGM neutral fractions (xjj i = 0.75, 0.50, 0.25), the bubble size increasing with decreasing neutral 
fraction. Note that ((tq/10^ yr) = 2.0) corresponds to Tq/rrec = 1.21 for xjj i = 0.5 and Tree xn , . The quasar is at the centre of each panel. The 
X-axis represents the LoS and the observer is on the right side of the panels. 
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0.5 
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0.79 
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-0.11 


0.29 
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46.3 


-0.03 


0.12 
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0.50 


49.3 


-0.04 


0.11 


4.0 


0.25 


52.8 


-0.06 


0.07 



Table 1. This tabulates the size (i?j_), and the anisotropy and shift param- 
eters (ri and s respectively) corresponding to the shapes shown in Figure 
2. 

not play a very important role if the quasar bubble is much larger 
than the typical size of the ionized patches outside the bubble, and 
we have ignored this contribution in our estimates in this paper. 
We have only considered the system noise ((Ai?)^)N which is the 
most dominant component in {{AE)^). The resulting SNR is in- 
versely proportional to the square root of the total observation time 
SNR <x tin . 

2.3 Anisotropic filter 

The matched filter technique considered in earlier works [Paper I; 
Datta et al. (2008); Datta, Bharadwaj & Choudhury (2009)] and 
briefly described above uses the signal expected from a spherical 
H II bubble [eq. (8)] as a template for the filter [eq. (11)]. The actual 
bubble is however anisotropic due to the FLTT, and we expect a 
mismatch between the template and the actual signal. We can avoid 
this if we use the apparent anisotropic shape as the template for the 
filter SfiU, v). This is done by numerically determining the bubble 
radius r as a function of (f) (Figure 1) by solving eq. (2) and (3). A 
section through the anisotropic bubble continues to be a circular 



disk, and eq. (8) still holds for the signal. The only difference is 
that the comoving radius Rb now varies along the LoS, and we 
have to use the calculated r(0) instead of a fixed Rt- An inspection 
of eq.s (1), (2) and (3) shows that the apparent shape is completely 
specified by three free parameters Nphs/C, xui /C and tq. It is 
therefore necessary to analyze the SNR for different values of these 
three parameters and determine the values of Npha/C, i /C and 
TQ for which the SNR peaks. 

Here we consider a targeted search around a known quasar 
whose infrared spectrum has been measured. It is possible to ex- 
trapolate the observed quasar's infrared spectrum blueward of the 
Ly-a line and estimate Nphs (example Mortlock et al. 2011). It 
is reasonable to assume that Nphs is known in a targeted search 
for bubble detection. Further, the analysis of numerical simulations 
indicate that the clumping factor has a value C = 30 at z = 8 
(Gnedin & Ostriker, 1997; Yu, 2005; Yu & Lu, 2005). With this as- 
sumption the apparent shape of the H II bubble is completely spec- 
ified by only two parameters the mean neutral fraction xh i and 
TQ the age of the quasar. Our entire analysis of bubble detection is 
in terms of these two parameters, i and tq. For a quasar with 
a constant luminosity {N ph s / sec~^) = 3, Figure 2 shows 
how the shape and size of the H II bubble vary with i and tq . 
The bubble grows as r-(r) oc {t /x-a i )^^^ in the early stages when 
T ^ Tree (cq. 2) and saturates at r(r) = Vs, which is indepen- 
dent of a;H I , when r/rrec ~ 1. Figure 2 and Table 1 show that the 
size of the bubble decreases with increasing xu i during the early 
stage of growth, and this effect becomes less noticeable at the later 
stage of growth. The bubble appears elongated along the LoS in the 
early stage, whereas it becomes compressed in the later stage. We 
see that the size anisotropy 77 and shift s all vary with tq and 

Xh I • 



3 ANALYTIC ESTIMATES 

In this section we present analytic estimates of bubble detection in 
a targeted search around a known high redshift quasar. For this pur- 
pose the H I outside the bubble is assumed to follow the dark matter. 
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Figure 3. This shows the analytic estimates of minimum observation time required for 3(j and 5(t (left and right panels respectively) detection of the H II bubble 
around the quasar ULASJ 11 20+0641 discovered by Mortlock et al. (2011). 



For the growth of the quasar bubble we ignore the effect of other 
ionizing sources and inhomogeneities in the H I distributions. The 
bubble's evolution is assumed to be exactly described by eq. (2) 
discussed earlier. For bubble detection we use the apparent shape 
obtained using eqs. (2) and (3) to calculate the filter Sf{lJ,v). 
The search has two parameters i and tq, and we expect a per- 
fect match when these have the same value as the bubble actually 
present in the data. 

Quasars typically have a luminosity of around Nphs ^ 
10^^ sec~^. The highest redshift quasar detected till date has a lu- 
minosity of Nph3 ~ 1.3 X lO^^sec^^ (Mortlock et al., 2011). 
Very little is known about the luminosity distribution of very high 
redshift quasars. The more luminous quasars are expected to pro- 
duce larger ionized bubbles, and the prospect of a 21 -cm detec- 
tion is also expected to be higher for a more luminous quasar In 
much of our analysis we have used two different quasar luminosi- 
ties Nphs = 3 X 10^^ sec~^ and 8 x 10^^ sec~^. The typical quasar 
age is expected to be in the range 10^ — 10* yrs (Haehnelt et al., 
1998; Haiman & Hui, 2001). The analysis of the proximity zones 
in the Ly-a spectra of high redshift quasars yield quasar ages which 
are larger than 1 - 3 x 10^ yrs (Worseck & Wisotzki, 2006; Lu & 
Yu, 201 1). We have considered quasar ages in the range lO'^' - 10* 
yrs in our analysis. 



3.1 Detectability 

We first consider the possibility of detecting the H II bubble around 
the quasar ULASJl 120-1-0641 (Mortlock et al., 201 1) using GMRT 
redshifted 21-cm observation. The size and shape of the H II bubble 
depend on xu i and tq , and we have determined the observa- 
tion time that will be required for different values of these pa- 
rameters. The left and right panels of Figure 3 show the obser- 
vation time that will be required for a 3cr and 5a detection re- 
spectively. We find that for 1,000 hr of observation a 3(7 detec- 
tion is possible in a considerable region of parameter space where 
x'H I > 0.4 and (tq /lO^ yr) > 1.5, and a 5(t detection is possible 
in a relatively smaller region of parameter space where xh i > 0.7 
and (rg/lO^yr) > 3.5 . A 3(7 detection is possible for a con- 
siderably large region of the parameter space (xn i > 0.2 and 



(rg/lO^ yr) > 1.0) with 4, 000 hr of observation, and we require 
approximately 8, 000 hr of observation for a 5a detection in this 
region of parameter space. 

The possibility of a 21-cm detection is quite favourable if the 
quasar luminosity is of the order of Nphs = 3 x 10^^ sec~^ or 
larger, and in much our analysis we use two different quasar lu- 
minosities Nphs = 3 X 10^^ sec^^ and 8 x 10^^ sec^^. Here we 
present estimates of the SNR for bubble detection considering 1000 
hr of GMRT observation around a known quasar. The right panels 
in Figure 4 show the peak SNR, which corresponds to the situa- 
tion when the filter and the signal are exactly matched. And the left 
panels in Figure 4 show the corresponding quasar bubble size (R± 
in Mpc). An ionized bubble grows as r(r) oc {rNphs/xn i ^^'^ in 
the early stages when r <^ Tree (eq. 2) and saturates at, 

= = 54Mpc(^^^^^)^'\ (15) 

for 

r > Tree = 3.3 X lO'' yr x^\ . (16) 

We see that both these stages are distinctly visible in the apparent 
size of the H II bubble shown in Figure 4. Note that the value of Vs 
where the growth saturates, is independent of xh i • 

The SNR from an H II bubble scales as x^ , Rl^^Nphs (Paper 

/ . \l/2 

I), which implies that SNR oc Ixm TNphs 1 when r <C Tree 

and SNR oc xu i -^V^^^ and independent of r when r > Tree ■ 

We observe that the SNR contours are roughly hyperbolas in 
the XHi — TQ plane, and the hyperbolas flatten out at large r. We 
mainly focus on the SNR= 3 contour, a 3a (or higher) detection is 
possible in the parameter region to the right of this contour We first 
consider the case where {Nphs /W^^ sec~^) = 3. We find that a 3(7 
detection is possible over a reasonably large region of the parameter 
space. For a high neutral fraction (0.6 < xui), it will be possible 
to detect an H II bubble even if it is small (R± ~ 20 Mpc) and in an 
early stage of its growth (rg/lO^yr > 0.5). However, it will not be 
possible to detect the bubble if the quasar's age {tq / lO^yr) is less 
than 0.5, even if the IGM is completely neutral outside the bubble. 
In contrast, a 3a detection is possible for a low neutral fraction only 



Constraining QSO and IGM Properties by Bubble Detection 7 



0.8 - 

0.6 - 

0.4 - 

0.2 i 








2 3 4 5 

tq/IO' yr 



2 3 4 5 

tq/IO' yr 



5s 



0.8 - 

0.6 - 

0.4 - 

0.2 I 








2 3 4 5 

tq/IO' yr 



2 3 4 5 

tq/IO' yr 



Figure 4. This shows analytic estimates of the bubble's comoving size iR± in Mpc){Left panel) and the SNR for the best matched filter (Right panel). Top 
and bottom panels correspond to quasar photon emission rates of {Nphs/W^'^ sec"'^) = 3 and 8 respectively. 



if the bubble is very large and in a later stage of its evolution. The 
lowest neutral fraction where a 3a detection is possible is a;H i ~ 
0.2. The quasar's age should be ~ 3.0 x 10^ yr (or more) for which 
the bubble radius is R± ~ 45 Mpc. For a quasar of this luminosity, 
a 3a detection is not possible if the neutral fraction is lower than 
0.2. We find that a 5a detection also is possible for a reasonably 
large region of the parameter space roughly covering xui > 0.3 
and (tq /lO'^yr) > 1. 

We next consider the case where (iVphs/10*^ sec~^) = 8. 
The scaling behaviour discussed earlier leads us to expect that R± 
and the SNR will increase by a factor of 1.38 and 1.63 respec- 
tively if {Nphs/W^''^ sec~^) is increased from 3 to 8. Note how- 
ever that, the scaling relations actually hold for the bubble seen 
from the rest frame of the quasar, and we do no expect them to 
be exactly valid for the bubble's apparent shape when the FLTT is 



taken into account. The effect of FLTT is expected to be larger for 
(A'phs/lO^^ sec~^) — 8 where the bubble is predicted to be larger 
in comparison to {Nphs/W"''' sec^^) = 3. 

We see that there is a increase in the region of parameter space 
where a 3a detection is possible for the more luminous quasar. The 
smallest age for which a 3a detection is possible is now reduced to 
(rg/lO^yr) = 0.3 compared to 0.5 when (iVphs/10^^ sec"^) = 
3. This however requires a completely neutral IGM. It is not possi- 
ble to detect the bubble if the quasar's age (Tg/lO^yr) is less than 
0.3, even if the IGM is completely neutral outside the bubble. At 
the other end, a 3a detection is possible if the neutral fraction is 
~ 0.1, provided the quasar's age is 2.5 x 10^ yr or more. For a 
quasar of this luminosity, it is not possible to detect bubbles if the 
neutral fraction is less than 0.1. 

There are indications that the neutral fraction could have a 
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Figure 5. This shows analytic estimates of the anisotropy parameter r) (Left panel) and the shift parameter s (Right panel) for the quasar's H II bubble. Top 
and bottom panels correspond to quasar photon emission rates of {Npf^^ /IQt''^ scc~^) = 3 and 8 respectively. 



value a-H i ~ 0.5 at 2 = 8 (Mitra et al., 2011). For this neutral 
fraction a Scr detection is possible for (rg/lO^yr) > 0.75 and 0.5 
for {Nphs/\0^'' sec~^) = 3 and 8 respectively. The corresponding 
H II bubble radii are 25 and 30 Mpc in these two cases respectively. 
We find that a 5(t detection is possible provided that the quasar's 
age exceeds (rg/lO^yr) > 2.0 and 1.0 respectively in these two 



We note that, the peak SNR does not increase very signifi- 
cantly when we use the anisotropic filter instead of the spherical 
filter used in earlier works. The anisotropic filter however has the 
advantage that it parametrises the bubble in terms of xh i and tq 
both of which are physically relevant and interesting quantities in 
their own right. We next consider the possibility of using matched 
filter technique to observationally determine xn i and tq. 



set 


I 


(yr) 


aO 


0.75 


1.0 


bO 


0.75 


1.8 


cO 


0.75 


2.8 


do 


0.50 


1.2 


eO 


0.50 


2.8 


fO 


0.25 


2.8 



Table 2. This tabulates the quasar and IGM parameters for which we have 
analytically considered parameter estimation. 



3.2 Parameter Estimation 

Once an H II bubble has been detected in 1, 000 hr of GMRT obser- 
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Figure 6. This shows the uncertainty (A SNR = 1) in the estimated parameters (xjj i and rg) for the parameter sets tabulated in Table 2. The three shaded 
regions correspond to 1, 000, 4, 000 and 9, 000 hr of observations respectively. The quasar photon emission rate is sec~^) = 3 in all the panels. 
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of using this to estimate the quasar's age and the IGM neutral frac- 
tion, TQ and xh I respectively. In most cases the 1, 000 hr obser- 
vation considered earlier is barely adequate for a detection, and we 
do not expect very significant parameter estimation from such an 
observation. Here we also consider deeper follow up observations 
of 4, 000 and 9, 000 hr specifically carried out for the purpose of 
parameter estimation once the H II bubble has been detected. 

Parameter estimation is based on the idea that tq and xh i 
are uniquely imprinted in the H II bubble's apparent size and shape 
(Figure 2). Here we use i?x (Figure 4) to quantify the bubble's ap- 
parent size, and 77 and s, introduced earlier, to quantify the bubble's 
apparent shape. Figure 5 shows the behaviour of rj and s for dif- 
ferent values of the parameters tq and xh i . Due to the FLTT the 
light reaching a present day observer was emitted earlier from the 
back side of the bubble in comparison to the light coming from the 
front side. For a bubble in the rapid stage of growth (r ^ Tree) the 
back side appears to have a much smaller radius compared to its 
front side, causing the bubble to appear elongated along the LoS. 
The centre of the apparent shape of the bubble also shifts towards 
the observer. Both these features are clearly seen in the parameter 
range (tq/ 10^ yr) < 1 (Figure 5) where the bubble appears elon- 
gated (77 > 0). We see that both 77 and s are extremely sensitive to 
TQ in this stage, however we do not see any significant xh i depen- 
dence. Though much of this region is outside the parameter range 
where a Scr detection is possible for 1000 hr of GMRT observation. 
For (iVphs/10^^ sec^^) — 3 and 8, there are small regions with 
high neutral fraction (xh i > 0.65) and (rg/lO^yr) > 0.5 and 
0.3 respectively where it may be possible to detect an elongated 
bubble. 

Most of the region ((tq/IO^ yr) > 1) of the i — tq pa- 
rameter space where the H II bubble is detectable, corresponds to 
the late stages of its growth where the bubble appears compressed. 
The bubble is no longer in the phase of rapid growth, and it appears 
mildly compressed (77 < 0) along the LoS. There is however still a 
considerable shift in the apparent centre of the bubble in this stage 
of growth. In this stage the anisotropy 77 is sensitive to both xh i 
and TQ . The shift s is mainly sensitive to tq and largely insensitive 
to Xh I (Figure 5). The bubble's radius (Figure 4) also is sensitive 
to both Xh I and tq in this stage of the bubble's growth. The pe- 
riod 0.6 — 4.0 X 10^ yr which is just after the rapid growth of the 
bubble appears to be the most suitable for parameter estimation as 
the apparent size and shape of the bubble are sensitive to both tq 
and Xh i • The bubble's apparent shape and size are only weakly 
sensitive to tq beyond (rg/lO^yr) > 4.0. However, the bubble's 
size remains sensitive to xh i , and it may be possible to constrain 
this parameter in the very late stage of growth. 

Ideally we expect the SNR to peak when the parameters (xh i , 
tq) of the filter exactly match those of the bubble that is actually 
present in the data. However the statistical fluctuations will cause 
the position of the peak to shift introducing an uncertainty in the 
parameter estimation. We use the criteria ASNR = 1 to estimate 
the uncertainty in the estimated values of the parameters xh i and 
TQ . Table 2 shows the sets of parameters for which we have consid- 
ered parameter estimation in this section. The parameter sets aO, bO 
and cO have a high neutral fraction (xh i = 0.75), dO and eO have 
Xh I = 0.5 and fO has a low neutral fraction (xh i = 0.25). Figure 
6 and 7 show the uncertainties in parameter estimation for all these 
cases assuming that (A'phs/lO^^ sec^^) — 3 and 8 respectively. 

There are indications that the neutral fraction could have a 
value xhi ~ 0.5 at z = 8 (Mitra et al., 2011). We have con- 
sidered two sets of parameter values (dO and eO) for this neutral 
fraction. Considering 1, 000 hr of observation for the set dO, where 



the bubble is in its early stage of growth ((rq/lO^yr) — 1.2), 
we find that it is possible to put both lower and upper limits on 
the quasar age (0.5 < (rg/lO^yr) < 3.5) for a quasar with 
(A'phs/lO^^ sec^^) — 3. In this case no limits on xh i can be 
placed with 1, 000 hr of observation. The limits on tq improves 
further and a lower limit can be placed on xh i when the ob- 
servation time is increased to 9, 000 hr. These limits now are 
1-0 < (rg/lO^yr) < 2.0 and xh , > 0.2. We see that there is 
a considerable improvement in parameter estimation if we have a 
brighter quasar with (iVphs/lO^^ sec^^) — 8. We now have the 
limits 0.5 < (rg/lO^r) < 2.5 and xh i > 0.1 for 1000 hr of 
observation, and 1.0 < (rg/10''yr) < 1.75 and xh i > 0.25 for 
9000 hr of observation. 

We next consider the set eO which corresponds to a later stage 
of growth ((rg/lO^ yr) — 2.8). We see that in this case the con- 
straints on tq and xh i are weaker compared to the set d where the 
quasar bubble is in an earlier stage of its growth. Here it is possi- 
ble to place only a lower limit on tq ((rg/10^ yr) > 1.5) for a 
quasar with (iVphs/ 10^^ sec~^) = 3 and 1, 000 hr of observation. 
These limits become 2.0 < (rg/lO^yr) < 5 and xh i > 0.15 
for 9, 000 hr of observation. In case we have a brighter quasar 
((A^phs/lO" sec^^) = 8), the limits are 1.75 < (rg/10^ yr) < 5 
and 2.25 < (tq /IQ' yr) < 3.75 for 1, 000 and 9, 000 hr of ob- 
servation respectively, and xh i can only be constrained ( 0.25 < 
Xh I ) with 9, 000 hr of observation. 

We next consider the possibility of a high neutral fraction 
(a;H I = 0.75) for which aO, bO and cO respectively correspond 
to progressively increasing quasar age. Comparing the results with 
those for xh i = 0.5, we find that in it is possible to place tighter 
constraints on tq and xh i if the neutral fraction is higher. We 
find that it is possible to place an upper and a lower limit on 
TQ , and a lower limit on xh i in all the cases that we have con- 
sidered.The constraints are tightest if the bubble is in the early 
stage of its growth, and we have 0.5 < (rg/lO^yr) < 2.25 
and Xh I > 0.1 for 1, 000 hr of observation for a quasar with 
(iVphs/lO^^sec^^) = 3 and (rg/lO^yr) = 1.0. The constraints 
are further improved if we have longer observations or a brighter 
quasar. 

We finally consider the set fO which corresponds to a low neu- 
tral fraction (xh i = 0.25). In this case the bubble is not large 
enough for a detection in the early stage of its growth (Figure 4), 
and we have considered (^g/lO^yr) = 2.8 which is in a later 
stage of its growth. We find that the prospects of constraining tq 
and Xh i are worse in comparison to the situation where we have a 
high neutral fraction. For a quasar with (Afphs/10*^ sec"^) = 3 it 
is not possible to constrain xh i even with 9, 000 hr of observation, 
and it is possible to place only a lower limits (rg /lO^ yr) > 1 and 
(rg/10^ yr) > 1.75 with 1, 000 and 9, 000 hr respectively. For a 
bright quasar ((A'pha/lO^^ sec~^) = 8), it is possible place con- 
straints 2 < (rg/10^ yr) < 4 and xh i < 0.75 with 9, 000 hr of 
observation. 

In summary, we find that the situation is most favourable for 
constraining tq and xh i if we can detect a very luminous quasar 
in the early phase of its growth and also in the early stage of reion- 
ization when the neutral fraction is high. 

As mentioned earlier in this section the anisotropy in the 
H II bubble's shape is determined by three parameters Nphs/C, 
Xh I /C and rg. We have assumed C = 30 for both the filter and 
the bubble. It is quite possible that the actual clumping factor is dif- 
ferent from C = 30. This will introduce a systematic uncertainty in 
parameter estimation. We consider two situations where the actual 
value of C for the bubble is different (C = 20 and 40) than the 
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Figure 8. This shows the estimated parameters in a situation where the value of C assumed for the filter is different from the actual value in the IGM. Top 
and bottom panels correspond to the IGM having C = 20 and 40 respectively. Left and right panels correspond to (tq/10^ yr) = 1.2 and 2.8 respectively. 
The filter assumes C = 30 in all cases. The bubble and the filter both have (A^phs/lO^^ sec~^) = 8. The square in each panel shows the estimates from 
the search while the circle shows the actual input parameters. The shaded regions show the uncertainty (A SNR = 1) in the estimated parameters for 1, 000, 
4000 and 9, 000 hr of observation. 



filter (C = 30) but the bubble and the filter both have same pho- 



ton emission rate {{Npha/^Q' 



— 8) in their models. Top 



and bottom panels of Figure 8 (corresponding to C = 20 and 40 
respectively) show the estimated values of parameters in these two 
different situations. We find that errors in C mainly affect the xn i , 
TQ is largely unaffected by this. In a situation where the C in the 
IGM of the bubble is less {C — 20) than what is assumed for the 
filter (C — 30), the filter severely underestimates xn i and this un- 
derestimation becomes more in the late stage ((rg/lO^ yr) = 2.8) 
of the growth than the early stage ((rg/lO^yr) = 1.2). On the 
other hand when the C in the bubble's IGM is more (C — 40) 
than that of the filter, the xh i is severely overestimated and this 
overestimation is more in the late stage ((tq/IO^ yr) = 2.8) of the 
bubble's growth. 

Similarly Nphs estimated from the spectrum of the quasar 
could be off from the actual value. It will also have similar effects 
as the errors in C in parameter estimations. 



4 SIMULATING THE IONIZATION MAP 

For our analytic estimates in the previous section we have consid- 
ered that the H I density is uniform in the IGM around the quasar 
and there is no other ionizing source in the IGM except the target 
quasar. It is expected that in reality the IGM H I density will not 
be uniform and the stellar sources in the vicinity of the quasar will 
also contribute in ionizing the IGM. To study the detectability and 



parameter estimation in a more realistic situation when all these ef- 
fects are taken into account we have simulated the ionization map 
around the quasar using a semi-numerical formalism. 

We have implemented the semi-numerical formalism pro- 
posed by Choudhury, Haehnelt & Regan (2009) for generating 
the ionization field at a given redshift. This uses an excursion-set 
formalism as introduced by Furlanetto, Zaldarriaga & Hernquist 
(2004). 

It is currently believed that stars residing in galaxies are the 
major source of photons to reionize the universe (Choudhury & 
Ferrara, 2006; Choudhury, 2009). The early stages of reioniza- 
tion are driven by stars within halos of mass ~ 10** h-^ Mp.. As 
the IGM becomes ionized, star formation within the smaller ha- 
los (M < ICp Mq) is inhibited because of radiative feed- 
back. Hence, it is sufficient to only include halos of mass M > 
lO'"* Mq to simulate the final stages of reionization. The small- 
est halo that is resolved in the simulation should be of mass M < 
10^ Mq. If at least 10 particles are required to constitute the 
smallest halo, then the particle mass should be < 10* h^'^ Mq. 
This decides the mass resolution of our simulation. 

We have generated the dark matter distribution at z = 8 
using a Particle Mesh Nbody code developed by Bharadwaj & 
Srikant (2004). The volume of the simulation is constrained by the 
16 Gigabytes of memory available in our computer. We perform 
our simulation in a periodic box of size 85.12 Mpc (comoving) 



with 1216"^ grid points and 608"^ 

Mpart = 7.275 X lO'^ h~'^ Mq. 



particles, with a mass resolution 
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Figure 9. This shows a simulated quasar H Ilbubble with sec~^) = 1.3 and xjj i = 0.5 at different stages of its growth, as seen from the 

quasar's rest frame. The square in each panel shows the location of the quasar. The expected spherical size is shown by the circle. Panels are labelled with the 
con'esponding quasar age and <p is the angle with the LoS as defined in Figure 1 and eq. (3). 



We identify halos within the simulation box using a standard 
Friend-of-Friend algorithm (Davis et al., 1985), with a fixed linking 
length fe = 0.2 (in units of mean inter particle distance) and mini- 
mum halo mass = lOMpart- We also compare the comoving num- 
ber density of halos per unit logarithmic mass interval dn/d In M 
with the theoretical mass function at z = 8 as predicted by Sheth 
& Tormen (2002) using the fitting function adopted from Jenkins 
et al. (2001). A good agreement is found over a wide mass range 
10^ < M/(/i"^M0) < 5 X 10". The lower mass limit corre- 
sponds to the smallest halo mass. 

The relation between the ionizing luminosity of a galaxy and 
its properties is not well known from observations. In the semi- 
numerical formalism adopted here, we assume that the ionizing lu- 
minosity from a galaxy is proportional to the mass of its halo. The 
number of ionizing photons contributed by a halo of mass M is 
given by. 



N-,iM) = N^on 



M 
niH 



(17) 



where niH is the mass of a hydrogen atom and Nion is a dimen- 
sionless constant. The value of Nion is tuned so as to achieve the 
value xn i desired in the simulation. In the semi-numerical formal- 
ism, a region is said to be ionized if the average number of pho- 
tons reaching there exceeds average neutral hydrogen density at 



that point. We have used a 256^ grid with a grid spacing of 0.3325 
Mpc for simulating the ionization maps. 

The semi-numerical formalism provides a snapshot of the ion- 
ized IGM at a fixed instant of time ts- We now briefly discuss 
how we have incorporated a quasar as an extra ionizing source in 
the simulation. The quasar is assumed to have ionizing luminos- 
ity Npha and age r at the instant ts- We have identified the most 
massive halo (M ~ 5 x 10^^ h"^ A-Iq) in the simulation as the 
quasar's location. We also shift the entire box with periodic bound- 
ary condition so as to bring the quasar into the centre of the field 
of view (FoV). In order to incorporate the quasar's H II bubble it 
is necessary to provide the simulation with a N-,^qso, the number 
of photons corresponding to the quasar at the instant ts. We have 
used. 



4 3 
N^,QSO = -Tv{nH)xHir (r) . 



(18) 



where r(r) is determined using eq. (2), i.e. the H Ilbubble would 
have a comoving radius r(r) if the hydrogen is uniformly dis- 
tributed and there were no ionizing sources other than the quasar. 

The quasar bubble actually produced in the simulation differs 
from this because of inhomogeneities in the hydrogen distribution 
and the presence of other ionizing sources. Further the density de- 
pendent non-uniform recombinations also change the bubble's size 
and shape in the cases where this effect is included. In general, we 
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Figure 10. This shows the apparent shape of the H II bubble around a quasar 
after the snapshots at different stages of its growth (Figure 9) have been 
staclced together. The square shows the location of the quasar The curve 
shows the expected apparent shape of the bubble. The observer is on the 
left of the figure. 



find the simulated H II bubbles are bigger than predicted by the eq. 
(2), mainly because of the contributions of other ionizing sources 
near the quasar. Figure 9 provides a visual impression of a suite 
of simulations with sec~^) = 1.3, xy_i = 0.5 and 

different values of r. We notice that there is a particularly large de- 
viation from the expected bubble size at the early stages of the bub- 
ble's growth. The actual bubble is much larger than that expected 
from the N-,^qso that we have assigned to the quasar. The massive 
halos, which are the only other ionizing sources besides the quasar, 
are preferentially clustered in the vicinity of the quasar which is lo- 
cated in the most massive halo in the simulation. The ionizing pho- 
tons from these halos outnumber those from the quasar in the early 
stages, and consequently the actual bubble is much larger than ex- 
pected. This discrepancy persists even in the later stages, however 
it is not so pronounced. A similar observation has been reported by 
Datta et al. (2012) in their radiative transfer simulations. 

A distant observer sees different parts of the H II bubble in 
different stages of its growth. We have assumed that the dark mat- 
ter distribution and the global xh i do not change within the look 
back time across the bubble. We have used eq. (2) and (3) to deter- 
mine r as a function of 0. We then choose sections each through a 
different member of the suite of simulations (Figure 9). Each sec- 
tion corresponds to a different value of r and t as determined by 
eq. (2) and (3). These sections are stacked together to produce the 
bubble's apparent shape as seen by the distant observer. Figure 10 
shows the apparent shape generated by stacking sections selected 
from the Figure 9. 

Our simulations spans a redshift interval Az = 0.28 along the 
LoS. We have made a simplifying assumption that the dark matter 
distribution and the variation in xn i across our simulation box can 
be ignored. Analytic estimates (for details see Appendix A and Re- 
sults of Paper II) show that the effects ignored here make a 5% 
or less contribution. Further, simulations (Datta et al., 2012) also 
reveal similar findings, justifying the assumptions made here. 



set 






tq/W 




(sec"-*-) 




(yr) 


A 


1.3 


0.50 


2.00 


B 


1.3 


0.50 


3.00 


C 


1.3 


0.75 


1.50 


D 


1.3 


0.75 


3.50 


E 


3.0 


0.50 


0.75 


F 


3.0 


0.75 


0.50 



|7 



Table 3. This tabulates the quasar and IGM parameters for which the 
H II bubbles have been simulated. 



The GMRT FoV at 151 MHz has a full width at half max- 
ima (FWHM) — 2.28° while our simulation box only subtends 
an angle 0.53° on the sky and 4.8 MHz along the LoS at redshift 
z = 8. Thus the simulation box will not be able to replicate the full 
GMRT FoV and can accommodate a bubble of maximum radius 
40 Mpc, which subtends ~ 0.5°. Maselli et al. (2007) predict 
the comoving radius of quasar generated H II regions to be ~ 45 
Mpc at 2 = 6.1 with xu i = 0.1. The size is expected to be less at 
z = 8. Our simulation box is thus large enough to host H II bubbles 
in the relevant size range at this redshift. The simulation adequately 
replicate the H I fluctuations in the vicinity of the bubble, however 
the H I fluctuations at large angular separation from the bubble's 
centre are missing. The H I fluctuations are some what underesti- 
mated as a consequence (see Datta et al. 2008 for more details). 

The simulated H I maps were converted to GMRT visibilities 
which were then used for matched filter bubble detection following 
the procedure detailed in Datta et al. (2008). 



5 RESULTS 

We have created six different realizations of the dark matter dis- 
tribution using the PM Nbody code. These were used to generate 
six independent realizations of the quasar bubble for each of the 
parameter sets given in Table 3. We have seen that it is possible to 
constrain different regions of the xui —tq parameter space by ob- 
serving quasars in different stages of growth. Ideally we would like 
to simulate bubbles in different stages of the growth for a variety of 
Nphs and xn i . However our choice of parameters are restricted by 
the condition that the bubble should be small enough (< 40 Mpc) 
to be contained within the simulation box, while it should also be 
large enough for a 3a detection (Figure 4). Figure 1 1 shows the 
simulated H II bubble for each set of the parameters listed in Table 
3. The images shown in the figure all correspond to the same real- 
ization, with uniform recombination. Three sets of simulations A, 
B and E have xn i = 0.5. The sets C, D and F have a higher neu- 
tral fraction (a;H i = 0.75). We are unable to simulate any situation 
with a low Xu i because the bubble size would have to be larger 
than our box for a detection to be possible. The sets A, B, C and D 
have a quasar luminosity (iVphs/lO^^ sec^^) = 1.3 which is the 
same as the quasar detected by Mortlock et al. (2011). This lumi- 
nosity is relatively low in comparison to the values that we have 
considered earlier. The quasar has to be observed at a late stage for 
a detection to be possible, and we expect a large bubble. The bubble 
centre will be shifted from the quasar for A and C where the quasar 
is somewhat younger with (tq / 10^ yr) = 2.0 and 1.5 respectively 
compared to B and D which have (rg/lO^ yr) = 3.0 and 3.5 re- 
spectively. The bubbles are also expected to be compressed in sets 
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Figure 11. This shows the simulated bubbles for different sets of quasar and IGM parameters as tabulated in Table 3 considering uniform recombination in the 
IGM. The square in each panel represents the location of the quasar. The curve shows the expected bubble shape. In all panels the observer is to the left side 
of the box. 



A and C, which is distinctly visible in these images. The bubble is 
expected to be nearly spheiical in sets B and D. 

We see that the simulated bubbles are in reasonable agreement 
with what is expected, though the simulated bubbles are nearly al- 
ways larger. This discrepancy arrises due to the extra ionizing pho- 
tons contributed by the halos located within the bubble. The num- 
ber of ionizing photons contributed by these halos is proportional 
to the ionization fraction (1 — a;H i )■ Consequently the discrepancy 
between the expected bubble size and that obtained in the simu- 
lations is larger for i = 0.5 in comparison to a;H i = 0.75. 
Further, the bubble's boundary is distorted because the bubble gets 
merged with the ionized regions produced by halos located in the 
periphery of the bubble. These distortions are particularly impor- 
tant when the halos located near the bubble's boundary contribute a 
significant fraction of the total ionizing photons within the bubble. 
We thus expect these distortions to be relatively important when 
the bubble is in the early stage of growth and we have a low neu- 
tral fraction, whereas we expect the distortions to be relatively less 
important when the bubble is in a later stage and in an IGM with a 
higher xm- 

The sets E and F have (iVphs/lO^'^ sec~^) — 3, and the 
quasar's ages are smaller {{tq/IQ^ yr) — 0.75 and 0.50) in com- 
parison to sets A,B,C and D discussed earlier We expect the bub- 
bles to be elongated along the LoS in these two sets, and this is 
clearly visible in Figure 1 1 . However, the bubble's shape is severely 
distorted in E where the neutral fraction is less compared to F. We 



expect such distortions to severely affect parameter estimation. We 
see that the distortion present in E also persists for F, however the 
relative contribution is smaller. We may expect the impact of this 
distortion on parameter estimation to be less severe in F as com- 
pared to E. 

Figure 12 shows images from simulations with inhomoge- 
neous recombination. We see that the bubbles are larger in com- 
parison to the situation with uniform recombination. It is possible 
to understand this by noting that in these models the dense regions 
in the IGM remain neutral through enhanced recombination, and 
only the low density regions are ionized. Consequently, a larger 
volume has to be ionized in order to achieve a particular value of 
the mass averaged neutral fraction xn i . Further, the contrast be- 
tween the bubble and the IGM outside is reduced, and the bubble's 
boundary is more distorted. 

For a fixed set of parameters, we expect the simulated bubble 
to differ from one realization to the next. Figure 13 shows two dif- 
ferent realizations of the simulated bubble for the parameter set A 
in addition to the realization shown in Figure 1 1 . We see that in all 
the realizations the bubble appears larger than expected. The distor- 
tions, we notice, can vary significantly from one realization to the 
next. In particular, the right panel of Figure 13 shows a realization 
where an extraneous ionization source is so aligned that it causes 
the ionized bubble to appear elongated along the LoS. Recollect 
that the bubble is expected to appear compressed in set A. Such 
distortions are a source of concern for parameter estimation. 
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Figure 12. Similar to Figure 1 1 considering non-uniform recombination in the IGM. 




Figure 13. Tliis shows two different realizations of the simulated H II bubbles having the same set of parameters (set A). The effect of the photon contribution 
from neighbouring halos is severe in the right panel as compared to the bubble in the left panel. The curve shows the expected bubble shape. In both panels 
the observer is to the left of the box. 
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Figure 14. This shows the SNR for 1, 000 hr of observation for the different sets of parameters in Table 3. The squai'es indicates the analytic estimates 
and the circles represent the SNR estimated from different realizations of the simulations. Left and right panels correspond to simulations with uniform and 
non-uniform recombinations respectively. 
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Figure 15. This shows the bubble size estimated from matched filter search on simulated H II bubbles for the different sets of parameters in Table 3. The 
squares indicate the analytic estimates of and the circles represent the same estimated from different realizations of the simulations. Left and light panels 
conespond to the simulations with uniform and non-uniform recombinations respectively. 



Distortions in the H II bubble and fluctuations in the 
H I distribution outside the bubble will affect the SNR for bubble 
detection and also the radius R± estimated from the best matched 
filter. Figure 14 shows a comparison of the expected SNR from our 
analytic estimates and the SNR obtained from the simulations. We 
note that the expected SNR has values 2.46 and 2.26 for the param- 
eters sets E and F, and a 3cr detection is not possible with 1, 000 hr 
of observation. We need 1, 500 and 1, 750 hr of observation for a 
3(T detection for E and F respectively. 

We find that the SNR is below the expected value in cases A 
and B, whereas it is comparable to the expected value in D and 
E, and it is larger than the expected value in C and F. We note 
that the reduction in SNR seen in A and B may be an artifact of 
the limited volume of our simulation where the bubble becomes 
comparable to the box size. The sets C and F where the SNR is 



larger than the expected value both correspond to an early stage of 
the growth in an highly neutral IGM (sh i = 0.75). We also find 
that there is a significant scatter in the SNR values amongst the 
different realizations. This scatter is relatively smaller for the sets 
which have a higher neutral fraction x^i =0.75 where we expect 
the halos to make a relatively smaller contribution to the ionizing 
photons. The scatter in the SNR remains the same even when we 
consider non-uniform recombination. However, the overall values 
of the SNR obtained from the different realizations are now slightly 
reduced with respect to uniform recombination. The fact that the 
SNR goes down if we introduce non-uniform recombination is a 
consequence of the lower contrast between the ionized bubble and 
the IGM outside. 

Figure 15 shows the bubble radius _Rx estimated from the best 
matched filter for the simulations compared with the analytic es- 
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Figure 16. This shows the location of the peak SNR for the different sets of parameters in Table 3. The squai'e in each panel indicates the actual value of xjj i 
and tq , the circles and stars correspond to the estimated values from different reahzations of the simulations with to uniform and non-uniform recombination 
respectively. The shaded regions coirespond to the analytic estimates of the uncertainty ( ASNR = 1) for 4, 000 and 9, 000 hr of observations. 



timates. As discussed earlier, we find that the simulated bubble is 
nearly always larger than what is expected from estimates. This dis- 
crepancy is particularly noticeable 15 — 20%) in the cases which 
have a lower neutral fraction xh i = 0.5 and where we expect the 
halos to make a relatively larger contribution. Further as noted ear- 
lier the discrepancy is larger in the simulations with non-uniform 
recombination in comparison to those with uniform recombination. 

We next use the simulated H II bubbles to analyze parameter 
estimation. Figure 16 shows the parameter values corresponding to 
the peak SNR for the different realizations of the simulations. We 
see that the values of tq and xh i recovered from the simulations 
are different from the actual age and neutral fraction. Reason for 
these deviations are the distortions in the bubble due to the other 
ionizing sources and the inhomogeneities in the IGM. Figure 16 
also shows analytic predictions of the uncertainties (ASNR = 1) 
in the estimated values of xh i and tq for the different sets (Table 
3) considering 4, 000 and 9, 000 hr of observation. We have not 
shown 1, 000 hr observation as bubble detection is not possible in 
this observation time for a few of the cases (E and F) that we have 
considered. We find that most of the simulations give parameter es- 
timates which lie within the region corresponding to ASNR = 1 
for 9, 000 hr of observation. We first discuss the sets E and F where 
the bubble is at a very early stage of its growth and is rather small. 
We require more than 1, 000 hr of observation for a 3(t detection. 
We find that in these two sets the neutral fraction is underestimated 
in almost all the realizations. This is a consequence of the fact that 
the actual bubble is almost always larger than expected (Figure 15) 
due to the presence other ionizing sources within the bubble. This 
effect is more pronounced when the IGM neutral fraction is low 



(a;H I = 0.5 in set E). This discrepancy in the bubble size only af- 
fects the estimated neutral fraction. The bubble in its early stage is 
very sensitive to tq through the anisotropy 77 and the shift parame- 
ter s (Figure 5) which are not severely affected by the overall size 
of the bubble. We find that in sets E and F the quasar age inferred 
from the simulations is quite close to the actual value. 

We next consider A, B, C and D which correspond to later 
stages of the bubble's growth. We find that the scatter in the val- 
ues of xh I and tq inferred from the simulations is considerably 
smaller for a high neutral fraction xh i = 0.75 (C and D) as com- 
pared to Xh I = 0.5 (A and B) where the halos make a relatively 
larger contribution to the ionizing photons. In sets A and B the es- 
timated Xh I is found to be uniformly scattered around the actual 
value. However, the value of tq is overestimated in the majority of 
the simulations. In set D we find that both xh i and tq are slightly 
overestimated in the simulations, while nearly all the simulated val- 
ues lie very close to the actual value in set C. 

In summary, our analysis is based on the assumption that the 
SNR of the estimator will peak when the parameters of the filter ex- 
actly match those of the bubble that is actually present in the data. 
The statistical uncertainty of the estimator will, however, cause the 
peak to shift introducing an uncertainty in parameter estimation. 
Further, distortions in the shape and size of the bubble due to other 
ionizing sources and the inhomogeneities in the IGM will also 
cause the peak SNR to shift, introducing further uncertainties in 
parameter estimation. In our work we have used the criteria ASNR 
= 1 to analytically estimate the statistical uncertainties. Further, we 
have used simulations to assess the effect of the distortions. We find 
that the statistical uncertainties are large in most of the situations 
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that we have considered, and hence these can be used to determine 
the uncertainties in the estimated parameter. We find that reliable 
parameter estimation is possible using H II bubbles for which a 3a 
detection is possible in 1, 000 hr of observation. Smaller bubbles, 
which require longer observations, do not provide very reliable es- 
timates of a;H I , they can however be used to obtain a reliable esti- 
mate of TQ . 

At the moment we do not have a rigorous justification for the 
criteria ASNR = 1 which we have adopted here. A more detailed 
statistical analysis involving extensive simulations with large num- 
ber of realizations of both the statistical fluctuations as well as the 
distortions are required. At the moment this is beyond the com- 
putational power available to us. However, 21 -cm observations of 
ionized bubbles hold the promise of allowing us to probe the age 
of the quasar and the neutral fraction of the IGM, and we plan to 
pursue such simulations in future. 
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